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The low temperature regime of non-Abelian two dimensional ferromag- 
nets is investigated. The method involves mapping such models into certain 
site-bond percolation processes and using ergodicity in a novel fashion. It is 
concluded that all ferromagnets possessing a continuous symmetry (Abelian 
or not) exhibit algebraic decay of correlations at sufficiently low tempera- 
tures. 
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In a recent letter [1] Seller and I proposed studying the phase structure of the 2D 
0{N) models by mapping them Into a correlated site-bond percolation problem. This 
approach was applied to certain discrete spin modes and to the 0{2) model, for which we 
rederlved the Proehllch and Spencer [2] result regarding the existence of a massless phase 
at sufficiently low temperatures 1/ f3. In this paper I report an extension of the percolation 
approach to 0{N) N > 3. It leads to the conclusion that a massless phase exists In all 
0{N) models. 

For completeness I will repeat the main points of Ref. [1] (see also Ref. [3] for a more 
complete discussion). With any 0{N) spin configuration one can associate an Islng spin 
configuration by dividing the sphere S{N — 1) Into two hemispheres and Introducing an 
Islng variable cr = ±1, which specifies In which hemisphere the spin points. In this manner 
the standard nearest neighbor action (s.n.n.a.) for the 0{N) model allows rewriting the 
partition function as 




ds\udspi ■ exp 



Spi ■ Spj) 



(1) 



Here u Is the unit vector chosen for specifying the hemispherical decomposition, = \s-u\ 
and Sp ■ u = 0. With respect to the Islng variables the action Is ferromagnetic, hence 
amenable to the Fortuln-Kasteleyn transformation [4]. This procedure associates to the 
Islng problem a correlated site-bond percolation process defined as follows: 
FKl-identify clusters of like-cr spins (H-clusters) 

FK2- within each H-cluster occupy bonds randomly with probability 1 — exp{—2Pij) 
(obtain FK-cluster) 

FK3-asslgn to every site within a given FK-cluster the same a value, obtained by 

choosing randomly + or - with probability 1/2. 
Here Pij is the space dependent Inverse temperature, which for the s.n.n.a. would be 
/9s||jS||j. Fortuln and Kasteleyn proved that the mean FK-cluster size (expected size of the 
cluster attached to the origin) equals the magnetic susceptibility of the Islng variable 



Xls 



^ E ^^-^^) • (2) 



In particular the latter diverges when the mean FK-cluster size diverges. 
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To apply the F-K procedure to the 0{N) models, Seller and I considered a modified 
model called 'cut' action: the Gibbs factor is s.n.n.a. only ii \si — Sj\ < e, < e < 2 and 

otherwise. We then formulated the following three conjectures: 
CI: The Mermin- Wagner theorem applies to the 'cut' model. 
C2: The 0{N) models ('cut' or not) are ergodic. 

C3: On a triangular lattice T a percolation process produced by a measure enjoying the 
symmetries of the lattice can contain at most one percolating cluster. 

1 refer the reader to Refs. [1] and [3] for a thorough discussion of the motivations behind 
these three conjectures and of the comparison of the 'cut' and the s.n.n.a. models. I will 
elaborate only on C2, which is central to the arguments presented in this paper. Imagine 
a very large lattice on which one has used the Monte Carlo procedure to simulate the 
0{N) model. If one has achieved thermalization, then this configuration is 'typical.' In 
the infinite volume limit a typical configuration has two important properties: 

PI: spacial averages equal ensemble averages (Birkoff's theorem) 
P2: the configuration is (statistically) invariant under additional Monte Carlo steps. 
I will briefly sketch the argument used in Ref. [1] to prove that the 'cut' 0(2) model 
must exhibit algebraic decay of its correlation functions for e sufficiently small. In Eq. (2) 
let (•) stand for expectation value measured with the full Gibbs measure. By PI and P2, 
Xi^ can be computed as a quenched expectation value provided the spins s^|| are assigned 
the values of a typical configuration. Since the Gibbs measure is invariant under lattice 
translations and (discrete) rotations, by Cl and C3 a typical configuration cannot contain 
a percolating H-cluster. An interesting theorem by Russo [5] states that if a translational 
invariant percolation process on a T lattice is such that neither clusters of the set E nor 
of its complement E percolate, then the mean cluster size of both E and E must diverge. 
Taking E to stand for a — +1 and E for a — —1 shows that the mean size of the H-clustcrs 
must diverge. (This statement is not surprising since at /? = and e = 2 the 0{N) model 
is equivalent to the Bernoulli site-percolation process with p = 1/2 and for the latter the 
critical density on a T lattice is indeed 1/2.) The FK-clusters are subclusters of the H- 
clusters obtained via rule FK2. In the 'cut' model, this rule must be amended. Indeed 
because of rule FK3, the constraint could be violated unless bonds are occupied at all sites 
having > d = e/2. Therefore, in a 'cut' 0{N) model, the FK-clusters must contain 
D-clusters defined by the condition > d. In Ref. [1] we showed that for the 'cut' 0{2) 
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model simple applications of Cl and C3 required that neither D-clusters nor D-clusters 
(s|| < d) can percolate and then, by Russo's theorem, both must have divergent mean size, 
QED. 

From the discussion presented thus far it follows that in any 'cut' 0{N) N > 1 
model on a T lattice, if neither clusters of D nor of D percolate, the mean FK-cluster 
size must diverge and hence correlations must decay algebraically. In fact in the 'cut' 
model D-clusters can never percolate. Indeed the set D consists of two disconnected 
pieces, both of which are contained in H-clusters and I have already argued that H-clusters 
cannot percolate. Therefore, the only question is whether 5-clusters could percolate for 
e sufficiently small? The reason for which a topological answer to this question exists in 
0(2) is that in that case the set D consists also of two disconnected pieces, which, for 
e < v^, cannot communicate. Obviously in 0{N)N > 3, l) is a connected set and a new 
strategy must be employed. In the sequel I will state three independent arguments, that in 
the 'cut' 0{N) model l)-clusters cannot percolate for e sufficiently small or (3 sufficiently 
large. Each argument requires a new conjecture and I will address their merits too. 
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Argument 1 

This is a proof by contradiction. For simplicity I will discuss the s.n.n.a. 0(3) 
model (e = 2) at large and choose u = z. I will take d small but independent of /3 - 
so that by FK2, when /3 is large, the bond occupation probability for sites in D goes to 
1. I will assume that a cluster of D percolates and show that that assumption suggests 
that a certain magnetic susceptibility (Eq. (4)) diverges. To that end I introduce spherical 
coordinates and rewrite the partition function as 

Z = I JJ^ / d6i d(pi j ■ exp < (3 ^ [cos 9i cos 9j + sin 9i sin 9j cos{(pi — cpj] > (3) 
VeA-^o Jo J [ ^. J 

Consider the following susceptibility 

^^ = JX\Y1 (COSI'/'X - ^y)) (4) 



By PI and P2 Xf could be measured by quenching the 9 variables to the values 9 they 
would take in a typical configuration. That is 



where {■)q means expectation value computed with the measure 



JJ / d(fi] exp 



/3 ^ sin 9i sin 9j ■ cos{(fi — (fj) 



(6) 



Since the quenched model is an 0{2) model (albeit with space dependent couplings), one 
can employ Ginibre's inequality [6] to bound Xf from below by the value it would take 
if in the measure (6) one replaced (3 sin 9i sin 9j by at all sites where a/ (3 sin 9i < c for 
some c > 0. Under the assumption that D percolates, by C3, these sites could not possibly 
percolate, but would form islands. The average size of these islands relative to the average 
distance between them would decrease with beta. Indeed by the Mermin- Wagner theorem, 
the probability of finding the spin at a site taking values in some subset of the sphere 
A of volume V{A) is equal to V{A)/A7r. (For P large, one can use perturbation theory 
to estimate the average size of these islands, which becomes actually independent of (3.) 
Thus the assumption that the equatorial strip D percolates implies that Xf is bounded 



from below by the susceptibility of an 0(2) model at large inverse temperature, but on a 
lattice having some small, randomly distributed holes. Although I am not aware of any 
rigorous result proving that, the following conjecture seems eminently reasonable. 

C4: Consider a T lattice and dilute bonds randomly with a probability smaller than the 
percolation probability for unoccupied bonds. Then there exists a Pkt < such that for 
any /3 > Pf-t the susceptibility diverges. 

Before motivating this conjecture, let me say that there is no reason to expect 
that if in the 0(3) model D percolated, the polar caps would be distributed as the holes 
produced by a Bernoulli process. Their actual distribution would be controlled by the full 
0(3) measure. However, if D percolated and especially if the model had a mass gap, by 
some central limit theorem, one would expect the polar caps to form islands and their 
distribution to be random at distances much larger than the correlation length. 

The intuition for C4 comes from the following rigorous results: 

a) Georgii [7] proved that if one randomly dilutes sites or bonds on a regular lattice 
with D > 2, then provided a remaining cluster percolates, there exists an inverse 
temperature f3c < oo such that for (3 > (3c there exists long range order (l.r.o.). 

b) De Massi et al. [8] , proved that under the same conditions as above, the Laplacian 
retains its continuous spectrum. 

In the language of the Coulomb gas, my conjecture is that if one introduces in the gas 
perfect conductors, randomly distributed, if the perfectly conducting regions do not per- 
colate, at sufficiently low temperatures, the Coulomb gas does not exhibit Debye screening 
(the introduction of the perfect conductors will only affect the dielectric constant). 

To conclude this argument, the contradiction is this: if one assumes that for the 
0(3) model D percolates and Xis finite, then clearly so is the s^-susceptibility (since 
Sz < 1). On the other hand C4 strongly suggests that the Sx — Sy susceptibility would 
diverge when P is large. This is a clear violation of 0(3) invariance, hence the assumption 
that D percolates must be false. Although not transparent, the topology of 0(3) is crucial 
for this argument. Indeed one may wonder if a similar reasoning could not be used to relate 
the 0(2) model to the Ising model and thus prove that the latter must exhibit l.r.o. at 
large /3, in violation of the Mcrmin- Wagner theorem? The answer is no, precisely because 
D is no longer a connected set and thus it could not possibly percolate. 
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Argument 2 

This is again a proof by contradiction. For simplicity I consider the 'cut' 0(3) 
model and choose u = z. I would like to argue that if the equatorial strip D percolated, 
then 0(3) invariance would be broken. 

Next let me consider the realistic case of a T lattice and an e small, yet e > 0. 
Suppose that in fact the equatorial strip D does percolate and hence its complement D 
forms islands. In the 'cut' model, the lines Sg = c> d will have to form closed loops, nested 
inside these islands. Consider now a c-tilted equator, namely the great circle passing thru 
Sz = c and Sx = 0. Since neither the hemisphere Sx > nor Sx < can percolate, any site 
of the lattice must be surrounded by an infinite sequence X{k) k e Z of concentric closed 
loops Sx = 0. (By the line = I mean a line on the dual lattice such that ■ s^. < 0; 
same type of qualifications apply to all other lines appearing in this discussion.) 0(3) 
invariance requires that the average number of intersections of the X lines with the c-tilted 
equators is independent of c. However if D percolates, then in any typical configuration 
there exists a /cq < oo such that any X{k) line with k > ko intersects the percolating 
cluster. That means that infinitely many X lines cross the c = tilted equator, while they 
may or may not cross the c-tilted equators with c > 0. In other words if D percolates, 
then one would expect the average number of crossings of the X lines with the c-tilted 
equators to decrease with c, in violation of 0(3) invariance. If on the other hand D does 
not percolate, then both D and D form rings and no a priori asymmetry in the average 
number of crossings of the X lines with the c-tilted equators exists. (An example where D 
percolates is the Richard model [9], which is a modified 0(3) model in which \sz\ < 1 — b 
for some 6 > 0, hence this model is only 0(2) invariant. The percolation approach used 
in Ref. [1] can be employed to prove rigorously that this model has to be massless for e 
sufficiently small - see Ref. [3]; Xy, diverges, yet Xis < oo-) 

In the discussion above I used the word 'expect' because one could say that even 
though if D percolates the regions with Sz > c are hidden inside regions of smaller Sg 
values, they are larger and thus restore 0(3) invariance. However 0(3) invariance requires 
that any typical configuration has the following two properties: 

Tl: The area is preserved. 

T2: The gradient is preserved. 
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Property Tl means that the density of sites where the spin points in some region A is 
proportional to the volume V{A). Property T2 says that if one selects two points on the 
sphere pi and p2, separated by a distance L, the average distance between sites where the 
spin points in the neighborhood of pi respectively p2 depends only on L (it is independent 
of which pi and p2 are chosen, provided they are at distance L). Obviously both properties 
are required by CI. 

C5: If in the 'cut' 0(3) model D percolated, then the typical configuration would violate 
Tl or T2 (or both) and, hence, 0(3) invariance. 

The motivation for C5 is this: if D percolated, then, as already argued, D would 
form islands - as opposed to rings, which are formed when neither D nor D percolates on 
a T lattice. The basic difference between a system forming islands and one forming rings 
is that islands are basically of finite size - the probability to find an island of diameter L 
decreases exponentially with L; on the contrary, if the system forms rings, there exists an 
infinite sequence of clusters surrounding each other and hence no exponential suppression 
of large clusters. Thus if the system forms islands the typical configuration will contain 
mostly mappings of a hemisphere over some finite region of T. It is easy to check that 
such maps cannot preserve both Tl and T2. No such difficulty exists if one considers rings 
- arbitrarily large regions of T. 

Argument 3 

As I have already noted, if D percolates, then D forms islands. Moreover, D consists 
of two disconnected pieces and Di- \Vhen e is sufficiently small, the volume of Dui 
V{Du) is much larger than that of D, V{D). On the other hand the area of the boundary 
of Du, S{Du) is half S{D). Is it reasonable to expect that under these circumstances, the 
mean cluster size of is finite while that of D infinite? The answer is provided by the 
following conjecture: 

C6: In the 'cut' 0{N) model, if two sets A and B have V(A)=V(B) and S{A) < S{B), 
then there exists eo{A, B) > such that for any e < eo, (A) > (B), where (•) represents 
the mean cluster size. 

The conjecture says that at given volume, the larger the surface of a set, the smaller 
its average cluster size. The reason for adding the qualifier that e < eo is that for e > the 
surface of the clusters of a set A need not consist of points on the surface of A. I believe that 
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this conjecture is intuitively clear. It can be proved in ID. In 2D it was verified numerically 
for 0(3) as follows: A was the Northern polar cap of area 47r/3, B the equatorial strip 
of the same area and e was such that the Northern and Southern polar caps could barely 
communicate. The data indicated that the mean cluster size of both A and B increased 
as L^~^ (L-linear size of the lattice) and that rjA < Vb- 

If C6 is true, it cannot be true that in the 'cut' 0{N) models the equatorial strip 
D percolates. Indeed if D percolates, its mean cluster size is divergent. By C6, for e 
sufficiently small, so is the mean cluster size oi D. By Russo's theorem, on a T lattice, 
that can occur only if neither D nor D percolates. QED. 
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Discussion 

The arguments presented above indicate that all 2D 0{N) models possess a mass- 
less phase. (This situation contradicts common wisdom. Evidence in favor of the latter 
is analyzed separately [10] and found wanting.) The arguments moreover suggest that 
although at large P extended topological defects - instantons - may exist in non-Abelian 
models, they are supressed entropically with respect to spin waves. This situation, already 
conjectured by the author in 1986 [11], suggests that for 

N >3 

the 2-point function may behave as 

(^o-^^)-<(^f-^+m^. (7) 

I have no basis at the present time to estimate a(/?) and b{P), nor whether rj depends on 
P in any given model. However, it could be that a and b are such that at intermediate 
distances the decay is exponential to a very good approximation (a similar effect governs 
the time evolution of a metastable state in nonrelativistic quantum mechanics [12]). 

Finally a word about perturbation theory. The fact that the 2D 0{N) models 
possess a massless phase for P suflBciently large does not imply that in 2D perturbation 
theory fails to produce the correct asymptotic expansion at fixed distances (as it does in 
ID for > 3). However if one defines the Callan-Symanzik /3- function by requiring that 
say (s(0) ■ s{x)) / (s(0) ■ s*(y)) is a renormalization group invariant for x, y 3> 1, then clearly 
an algebraic decay for fi > Pkt{N) implies that the Callan-Symanzik /^-function could be 
chosen to be vanishing. If my conjecture about Eq. (7) proved to be correct, one could also 
define the /^-function as dp/dln{m), in which case one may find the famous asymptotic 
freedom answer. However I find it hard to believe that if that were the case, the continuum 
limit constructed by letting — > oo would not contain (coupled) massless excitations (of 
course a continuum limit could also be constructed for any oo > /? > f3kt{N) - that field 
theory would be a massless theory) . 

Many of the ideas expressed in this paper stem from my long time collaboration 
with Erhard Seller. I am also gratefulfor the hospitality extended to me by the Max Planck 
Institut fur Pysik und Astrophysik - Munich. 
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